Double superconducting transition and low-temperature specific heat study of 

Pr(0 Sl ,,Ru :c ) 4 Sb 12 
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Fhe double superconducting transition in PrOs4Sbi2, first observed in specific heat C(T) mea- 
surements of single crystal samples, is the subject of an in-depth study. The double superconducting 
transitions of a batch of single crystals were measured before and after they were annealed for 5 days 
at 500°C, with no observed change. C(T) measurements near T c for PrOs4Sbi2 in several magnetic 
fields are also presented, detailing the evolution of the double transition up to 1 T. Samples of 
Pr(Osi_ I Ru a; )4Sbi2 with 0.01 < x < 0.04 also appear to display a double superconducting transi- 
tion in specific heat. In addition, samples with the smallest Ru concentration measured (x = 0.01) 
may even display a different type of superconductivity than pure PrOs4Sbi2 (x — 0). 

PACS numbers: 65.40.Ba, 71.27. +a, 74.25.Bt, 74.62.Dh 



I. INTRODUCTION 

The compound PrOs4Sbi 2 has attracted a great deal 
of attention in recent years as the first Pr-based heavy 
fermion superconductor, with a superconducting transi- 
tion temperature T c = 1.85 K and an electron effective 
mass m* « 50m e i£ Detailed measurements show the 
strongest support for a nonmagnetic T\ singlet ground 
state in a tetrahedral crystalline electric field (CEF), sep- 
arated from a T$ (in the cubic notation) triplet first ex- 
cited state by approximately 10 Ki 3 i 4 i 5 Quadrupolar or- 
dering arises below 1.5 K and above 4.5 T, related to 
the crossing of the Zeeman-split CEF energy levelsi&Sil 
The superconducting state appears to be highly uncon- 
ventional, exhibiting time-reversal symmetry breaking, 8 
multiple superconducting phasesf2il£ and double super- 
conducting transitionsi^iiiiiSiiii^ 

Upon substitution of Ru to form P^Osi-^Ru^Sb^, 
the transition temperature decreases to a minimum of 
T c « 0.75 K at x = 0.6, whereupon it increases to ~ 1.0 
K for PrRu4Sbi2. This behavior was observed in mea- 
surements of electrical resistivity, magnetic susceptibil- 
ity, and specific heat>i£*i& All three measurement tech- 
niques indicated that the CEF energy level splitting be- 
tween the ground state and first excited state increased 
monotonically from ~ 10 K for PrOs4Sbi2 to ~ 70 K for 
PrRu4Sbi2. The C(T) experiments also revealed that the 
electronic specific heat coefficient 7 decreased from ~ 500 
mJ/mol K 2 for x — to a plateau of 7 ~ 100 mJ/mol 
K 2 for x > 0.6. Upon subtracting the normal state 7 
and Schottky anomaly due to CEF effects, we were able 
to investigate the structure of C(T) below T c . We per- 
formed phenomenological fits which suggested that the 
PrOs4Sbi2 data below T c were more in agreement with 
power law behavior, indicating nodes in the energy gap, 
while the samples with 0.05 < x < 0.2 suggested expo- 
nential behavior, consistent with an isotropic energy gap. 
In addition, no obvious double superconducting transi- 
tion was observed for x — 0.05, the lowest substituted 
concentration measured in the previous study. 



In this paper, we present the results of further probes 
into the nature of the double superconducting transi- 
tion in PrOs 4 Sbi2. We observe the effect of magnetic 
fields, annealing, and Ru substitution on the double tran- 
sition through measurements of C(T). These results 
are compared to the double superconducting transitions 
observed in the heavy fermion systems U(Pti_ a; Pd 2: )3 
and Ui-zTh^Be]^. The low-temperature specific heat 
of Pr(Osi_ 2 ;Ru x )4Sbi2 is investigated in order to deter- 
mine a possible concentration at which the superconduct- 
ing state changes from one possessing nodes in the energy 
gap to one that is more isotropic. 



II. EXPERIMENTAL DETAILS 

The single crystals of Pr(Osi_ 2; Ru l: )4Sbi2 studied in 
this work were grown using an Sb flux method, and the 
crystal structure was confirmed via x-ray diffraction mea- 
surements, as described previously^ The specific heat 
C(T) measurements were performed in the same cryo- 
stat as for previous experiments^ using a semi-adiabatic 
method at temperatures between 0.6 K and 20 K. Mag- 
netic fields were applied using a 5 T superconducting 
magnet attached to the outside of the vacuum can. The 
experiments in a magnetic field were performed on the 
same individual single crystal measured by Cichorek et 
aLjifi with a mass of 8.17 mg. The field direction was 
along a principal crystal axis. The other measurements 
were made on collections of single crystals, with total 
masses between 20 and 30 mg. 

The annealing was carried out by wrapping the crystals 
in tantalum foil, which was sealed inside a quartz tube 
with a piece of zirconium foil, under 150 Torr UHP Ar. 
The quartz tube was placed in a 500 °C furnace for 5 
days, and then quenched to room temperature. The mass 
of the sample was the same after annealing within the 
systematic error of the weighing scales. 



2 



III. RESULTS AND DISCUSSION 

The only other stoichiometric heavy fermion supercon- 
ductor that displays a double superconducting transition 
is UPt 3 iA£ii£For single crystals of this compound, anneal- 
ing can be an important factor in observing the double 
transition in specific heatiiS The results of an anneal- 
ing study on a batch of PrOs4Sbi2 crystals are shown 
in Fig. 2] The two data sets fall neatly on top of each 
other, which argues against any kind of inhomogeneous 
mass redistribution. Unlike UPt3, however, the double 
transition in PrOs4Sbi2 shows no indication of becoming 
more resolved due to annealing. It is still important to 
note that the double transition does not degrade with 
annealing, as might be expected if this feature were not 
an intrinsic effect of PrOs4Sbi2. 

Figure|21displays the evolution of the double supercon- 
ducting transition with magnetic fields for PrOs4Sbi2. 
The double transition is clearly visible in the zero field 
data, while by 1 T the transition has been suppressed be- 
low 1 K. As the magnetic field is increased, the first tran- 
sition appears to lose its sharpness, until by 0.75 T it is 
either nonexistent, or indistinguishable within the strong 
downturn of the Schottky anomaly. The inset to Fig. [21 
shows the data after subtracting the 1 T curve, which fur- 
ther emphasizes the initial jumps. This evolution of the 
transition with magnetic field has been reported by sev- 
eral other groups, with mostly different results. Vollmer 
et al. found that the double transition disappears by 0.4 
Tji£ while the double transition observed by Grube et 
al. persists up to 0.6 TjSi similar to this work. The data 
reported by Measson et al. are noticeably different, as 
two clear transitions were visible up to ~ 2 T. 21 All re- 
sults, except those of Vollmer et al., were from experi- 
ments performed on one single crystal of PrOs4Sbi2; the 
experiments of Measson et al. relied on an ac calorimetry 
method instead of the semi-adiabatic method utilized by 
Vollmer et al., Grube et al., and this work. In UPt3, 
applying a magnetic field and lowering the temperature 
causes the two transitions to merge near ~ 0.5 T and 
~ 0.4 K, depending on the field orientation, giving rise 
to an additional phase transition within the supercon- 
ducting state^ 2 ^ While several experiments provide ev- 
idence for multiple phases in the superconducting state of 
PrOs4Sbi2, none of these phases can as yet be associated 
with the double superconducting transitionAifi 

Specific heat divided by temperature C{T)/T data for 
Pr(Osi_ x Ru x )4Sbi2 with < x < 0.05 are displayed in 
Fig. [3J The normal state data below 10 K were fitted 
with an equation that includes electronic, lattice, and 
CEF contributions, with the results listed in Tabled The 
quantity r is a scaling factor for the Schottky anomaly, 
as discussed previously^ This suppression of the Schot- 
tky anomaly could result from an energy dispersion due 
to Pr-Pr interactions or hybridization between the Pr f- 
electrons and ligand states, which results in entropy being 
transferred to the conduction electrons, thereby increas- 
ing the electronic specific heat coefficient 7. The errors 



for the parameters were determined by varying the Debye 
temperature 0d by ±10 K from its best fit value, as de- 
scribed previously^ Only results from a Ti ground state 
fit are presented here; we were able to accurately use the 
simpler cubic CEF equations since these measurements 
were only performed in zero field. As can be seen from 
Table [U the splitting between the ground state and first 
excited state, 6, increases monotonically with increasing 
x, as expected. The electronic specific heat coefficient 7, 
on the other hand, displays a maximum at x = 0.02. The 
Debye temperature 0d is also enhanced for x — 0.01 and 
x = 0.02 compared to the neighboring concentrations. 

The superconducting transitions in Fig. 0] are plotted 
as C e i(T)/T, where C e \ is the electronic specific heat after 
the lattice and Schottky terms were subtracted. Struc- 
ture in the superconducting transition is clearly apparent 
for x = 0.01 and x — 0.02, and also appears to be present 
for x = 0.04. In light of these data, the rounded structure 
of the x = 0.05 transition published previously^ might 
also hide a double transition, but it is too broad to tell 
for sure. The transitions for < x < 0.04 were analyzed 
using an entropy conserving construction for both transi- 
tions. The upper and lower transitions and their respec- 
tive specific heat jumps are presented in Table ITU It is in- 
teresting to note that the lower transition has the largest 
jump for PrOs4Sbi2, while for 0.01 < x < 0.04, the upper 
transition displays a larger jump. This is shown graph- 
ically in Fig. [3 which allows us to speculate that the 
upper transition jump and the overall transition jump 
may converge near x = 0.05. Also listed in Table UJ] 
is the ratio AC/jT c for the full transition. The value 
for x = of 1.76 is enhanced over the BCS value of 
1.43, suggesting strong-coupling superconductivity. The 
results for the upper and lower transitions are not tab- 
ulated, although they can be easily calculated from the 
other supplied data. The upper and lower transitions 
have AC/jT c values of 0.86 and 1.12, respectively. As x 
increases, AC/-fT c drops sharply, due to the strong in- 
crease in 7 and the smaller decrease in T c , to values well 
below 1.43, even for the full transitions. 

The evolution of the double transition in 
Pr(Osi_ 2; Ru a ;)4Sbi2 is qualitatively similar to that 
seen in U(Pti^. 2; Pd a ;)3, which displays a double transi- 
tion down to x = 0.003, but with no superconductivity 
observed at all for x > 0.004 above 1 K»Si The lower 
transition is very small in U(Pti_ 3; Pd 2 ;)3 for x = 0.003, 
but still may be present for higher concentrations at 
temperatures lower than what was measured. Never- 
theless, the specific heat jump at the lower transition 
in both series of compounds decreases in magnitude 
faster than the jump for the upper transition. An 
interesting discrepancy between the two substituted 
series is that in U(Pti_ a; Pd a: )3, the difference between 
the two superconducting transitions, AT C , increases 
with increasing x, while in Pr(Osi_ x Ru x )4Sbi2 the 
separation between the transitions remains nearly 
constant. Unfortunately, the smaller amount of research 
on PrOs4Sbi2 makes it difficult to determine which simi- 
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larities and differences between P^Osi-^Ru^Sb]^ and 
U(Pti_ a: Pd a ;)3 are meaningful. It should be noted that a 
double superconducting transition has been reported for 
Pri_ a; La a ;Os4Sbi2 with x — 0.02 and x = 0.05, and the 
lower transition also appears to be the one suppressed 
with La substitution^ 

Figure [fj] displays fits to the data below T c for both 
power law (appropriate for energy gaps with nodes) and 
exponential (appropriate for isotropic energy gaps) func- 
tions, a continuation of the fits performed previously^ 
The new fits for x = 0.01,0.02, and 0.04 were all made 
between ~ 0.6 K (the base temperature of the measure- 
ment) and 1.15 K. Neither fit is better than the other 
within the fit range, but the phenomenological extrap- 
olation of the fits to higher temperatures suggests that 
there may be a fundamental difference between the su- 
perconductivity in PrOs4Sbi2 and the Ru-doped mate- 
rials, even with only 1% Ru. A specific heat study of 
Ui-zTLj-Beia deep in the superconducting state revealed 
a similar evolution from power-law to BCS-likc behavior 
with increasing x, separated at x w 0.02.— The param- 
eters from fits to the Pr(Osi_ 2 ,Ru a ;)4Sbi2 superconduct- 
ing state data are listed in Table ITU A rough estimate 
for the coupling strength can be obtained by calculating 
2A e /T c (with A e in units of K), which is approximately 
4.5 for the values of x studied in this work. This en- 
hancement over the BCS estimate of 3.52 is consistent 
with other measurements on PrOs4Sbi2 which indicate a 
strong-coupling energy gap^i2& 

A double superconducting transition is one of the clear- 
est indications of unconventional superconductivity, and 
possibly a multi-component superconducting order pa- 
rameter. Studies of UPt3 strongly indicate that the lower 
transition in this compound arises from the onset of a 
new order parameter (presumably antiferromagnetic)^ 
the double transition in PrOs4Sbi2 is not, as yet, as 
well studied or understood. Time reversal symme- 
try breaking (TRSB) has been observed in PrOs4Sbi2 
through measurements of muon spin rotation (/j,SR) in 
zero field, which exhibit a spontaneous magnetic mo- 
ment arising in the superconducting stated It is diffi- 
cult to tell from the published data whether the moment 
arises due to the upper or the lower superconducting 
transition. In UPt3, a spontaneous moment in /iSR was 
reported below the lower superconducting transition, 29 
but this result has not proven to be reproducible via 



/j,SR measurements ^2i2i However, ultrasound, vortex, 
and flux flow measurements have revealed features be- 
low the lower superconducting transition associated with 
additional superconducting order parameters and, possi- 
bly, with TRSBii2i22i2ii2i A magnetic moment has been 
observed to arise in /xSR measurements below the lower 
superconducting transition of Ui-^Th^Beis with 0.02 < 
x < 0.04, also suggesting TRSB. 35 Measurements of vor- 
tex motion in Ui_ x Th x Bei3 with x = 0.0275 below the 
lower superconducting transition indicate a sharp drop in 
vortex creep rates, very similar to that seen in UPt3^& 
It therefore would be very useful to perform these mea- 
surements on low x samples of P^Osi-^Ruz^Sb^ in 
order to determine if the TRSB and unconventional su- 
perconductivity truly track the double superconducting 
transition, as these present measurements suggest. 



IV. SUMMARY 

In summary, we have investigated the double super- 
conducting transition in PrOs4Sbi2 through specific heat 
measurements in more detail than previously reported. 
Annealing the crystals did not appear to affect the fea- 
tures in the specific heat near T c . In a magnetic field, 
the transitions rapidly broaden in temperature and are 
suppressed by 1 T. Samples of Pr(Osi_ 2; Ru. r )4Sbi2 with 
smaller Ru concentrations than previously measured dis- 
played double superconducting transitions in C(T) up 
to x = 0.04. Power law and exponential fits to C(T) 
below T c also suggest that the superconducting state 
for x = is different than that for x > 0.01. The 
behavior of the double superconducting transition for 
Pr(Osi_ 2; Ru 2 ;)4Sbi2 appears to only share qualitative 
characteristics with the superconducting heavy fermion 
systems U(Pti_ a; Pd a; )3 and Ui-^Th^Beis. 
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FIG. 1: Specific heat divided by temperature C/T versus 
T below 10 K for an unannealed set of single crystals of 
PrOs4Sbi2 (open circles) and the same sample annealed for 
5 days at 500 °C (solid line). Inset: close-up of the dou- 
ble superconducting transition. There is almost no difference 
between the two data sets. 




FIG. 2: Specific heat divided by temperature C/T versus T 
between 1 and 2 K in magnetic fields up to 1 T for an in- 
dividual single crystal sample of PrOs4Sbi2. The data for 
fields below 1 T have been truncated below the supercon- 
ducting transitions for clarity. Inset: AC/T versus T, where 
AC = C(H, T) - C(H = IT, T). 
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TABLE I: Normal state physical properties of samples of 
Pr(Osi_ a: Ru :l ;)4Sbi2, determined from specific heat data. The 
parameter S is the splitting between the ground state and the 
first excited state in the Schottky anomaly, r is the scaling 
factor for the Schottky anomaly (as described in the text), 7 
is the estimated electronic specific heat coefficient, and 6d 
is the estimated Debye temperature. The errors in the pa- 
rameters were determined by allowing Od to vary by ±10 K 
within the fits (see text and RefJ- for details). 
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TABLE II: Superconducting state physical properties of sam- 
ples of Pr(Osi_2,Ru^)4Sbi2, determined from specific heat 
data. T c is the superconducting transition temperature, AC 
is the jump in C(T) at T c , n is the exponent of the power-law 
fit below T c , 7^ is the electronic specific heat coefficient in the 
superconducting state (7J = for the power law fits), and 
A c is the parameter in the exponential fit below T c that is 
proportional to the energy gap. The errors in the parameters 
were determined by allowing On to vary by ±10 K within the 
fits (see text and Ref^- for details) , with the exception of T c 
and AC/T C (the errors for AC/T C are represented graphically 
in Fig. . 
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FIG. 3: Specific heat divided by temperature C/T versus T 
below 10 K for Pr(Osi_ I Ru I ) 4 Sbi2 with < x < 0.05. In- 
set: superconducting transition temperature T c (filled circles, 
left axis) and the electronic specific heat coefficient 7 (filled 
squares, right axis) as a function of Ru concentration x. The 
vertical lines for T c delineate the upper and lower transitions 
for those concentrations which display double transitions. 




FIG. 4: Comparison of double superconducting transitions 
in Cci/T, after lattice and Schottky anomaly terms corre- 
sponding to a Ti ground state have been subtracted, for 
Pr(Osi_ a; Ru a; )4Sbi2 with (a) x = 0, (b) x = 0.01, (c) 
x = 0.02, and (d) x — 0.04. The transitions were approxi- 
mated with an entropy conserving construction. 
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FIG. 5: Specific heat jump at T c , AC/T C , for PrOs 4 Sbi 2 as 
a function of Ru concentration x. The dashed lines represent 
speculations on the continued evolution of AC/T C . 




FIG. 6: Comparison of exponential (dashed line) and power 
law (solid line) fits below T c , after subtracting lattice 
and Schottky anomaly terms for a Ti ground state, for 
Pr(Osi_o : Rua ; )4Sbi2 with (a) x = 0, (b) x = 0.01, (c) 
x = 0.02, and (d) x = 0.04. The x = fits only extend 
up to ~ 1.2 K, while the x = 0.01, x = 0.02. and x = 0.04 
fits only extend up to ~ 1.15 K, as described in the text. 



